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We reexamine the question of whether the geometrical ground state of neutral and ionized clusters 
are identical. Using a well defined criterion for being 'identical' together extensive sampling methods 
on a potential energy surface calculated by density functional theory, we show that the ground 
states are in general different. This behavior is to be expected whenever there are meta- stable 
configurations which are close in energy to the ground state, but it disagrees with previous studies. 

PACS numbers: 



I. INTRODUCTION 



II. METHODOLOGY 



Since experimental mass selection methods require ion- 
ized systems, the majority of experimental information 
on clusters was obtained for ionized clusters. On the 
other hand, neutral systems are of greater practical in- 
terest and the majority of theoretical works are done 
on neutral systems. The relation between the proper- 
ties of neutral and ionized clusters is therefore an impor- 
tant one. The basic property which determines all other 
properties is the structure. Finding the global minimum 
structure of a cluster is a complex global geometry opti- 
mization problem on a high dimensional potential energy 
landscape [l| with a huge number of local minima. In or- 
der to make accurate structural predictions, the potential 
energy surface should be calculated within density func- 
tional theory. Doing exhaustive unbiased searches for the 
global minimum at the density functional level has only 
recently become possible through the combined improve- 
ments in global optimization algorithms and computer 
performance. 

One basic question concerning the relation between 
neutral and ionized clusters is whether they have the 
same basic structure. Evidently adding or removing one 
electron will change the the exact bond lengths and an- 
gles but one might suspect that the structures remains 
nevertheless very similar. The relation between the struc- 
ture of neutral and ionized clusters has been investigated 
in numerous previous publications for the same silicon 
and magnesium clusters that we have reexamined. The 
conclusion, in all the publications we are aware of, is that 
in general the structures of the neutral and cation clus- 
ters are more or less identical, but the criteria for being 
'identical' are not always explicitly given. We introduce 
a well defined criterion for being identical. Two minima 
are identical or more precisely 'related', if the equilibrium 
structure of the ionized system lies within the catchment 
basin of the neutral system and vice versa. Applying this 
criterion on an extensive database of accurately relaxed 
geometries, we arrive at the opposite conclusion. 



The global and local minima presented here are ob- 
tained within Density functional theory using the'Big 
DFT' wavelet code 2 which was coupled to the 'minima 
hopping' Q global optimization algorithm. The local 
spin density approximation (LDA) is used together with 
HGH type pseudo potentials |4[ for the calculation of the 
potential energy surface. The size of the wavelet basis 
set was chosen such that the energies were converged to 
within better than 10 -4 Hartree with respect to the infi- 
nite size basis set. A combination of conjugate gradient 
and BFGS methods [E[ was used for the local geometry 
optimizations and they were stopped when the numerical 
noise in the forces was about 20 percent of the total force. 
This happened usually when the largest force acting on 
any atom was less than 2 x 10 -5 Hartree/Bohr. Saddle 
points were found by a modified version of the 'A spline 
for your saddle' method [6[. 

In contrast to plane wave basis sets, free boundary 
conditions for charged systems are not problematic with 
a wavelet basis set. In plane wave program a neutralizing 
background charge is needed, since a periodic system can 
not have a charged unit cell. In a wavelet basis set the 
integral equation for the potential V 
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can be solved directly for the electronic charge density 
p with a monopole and the electrostatic potential can 
therefore be calculated very accurately for charged sys- 
tems 0. 

For all the clusters we have carried out separate global 
optimization runs for neutral and ionized system. Since 
anions with weakly bound additional electrons are less 
accurately described by density functional theory than 
cations, we considered only cations in addition to the 
neutral system. For small clusters ( less than 10 atoms 
for silicon and less than 20 atoms for magnesium) the ma- 
jority of low energy local minima can be obtained. That 
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this condition is fulfilled can be deduced in the minima 
hopping algorithm from a strong increase in the kinetic 
energy of the molecular dynamics trajectories. For larger 
clusters this explosion of the kinetic energy [H can not 
be observed for any reasonable short simulation time. 
In case of medium sized clusters we calculated always 
at least 100 low energy local minima structures and we 
did various empirical checks to convince ourselves that 
the global minimum was found. We checked for instance 
always that the lowest energy structures found for the 
cation system did not relax upon addition of an electron 
into a structure that was lower in energy than the puta- 
tive global minimum found for the neutral system. 

Using this approach we investigate whether the global 
minimum structures of neutral and positively charged 
clusters are related. We will use the following two cri- 
teria as the definition for two structures of a neutral and 
ionized system to be "related" 

• The equilibrium structure i of the cation will re- 
lax into the equilibrium structure j of the neutral 
cluster when an electron is added. 

• The equilibrium structure j of the neutral cluster 
will relax into the equilibrium structure i of the 
cation when an electron is removed. 

By relaxations we mean local geometry optimization with 
a sufficiently small step size, which will make it very un- 
likely that the local geometry optimization jumps out of 
the catchment basin within which the local geometry op- 
timization was started. The structures of the neutral and 
ionized system are thus considered to be related, if there 
is a one-to-one mapping between the global minima struc- 
tures upon addition and removal of an electron. This 
definition of two structures being related is motivated by 
the fact that the removal or addition of an electron in 
an experiment is quasi instantaneous on the time scale of 
the motion of the heavy nuclei. A cluster will therefore 
relax experimentally into the minimum of the catchment 
basin in which it finds itself after the addition or removal 
of an electron. In order to see whether our definition 
is fulfilled or not, we have introduced mapping charts 
that show which local minimum of the neutral system 
relaxes into which local minimum of the ionized system 
and vice versa. We consider the global minima structures 
of the neutral and ionized cluster to be identical if the 
two global minima structures are related according to the 
above definition. 

In order to detect the degree of similarity between two 
structures with N at atoms and atomic coordinates R a 
and R b respectively we have also calculated the configu- 
rational distance D 
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permutated in the search for the smallest possible D. It 
turns out that structures, that are related according to 
our definition, usually have also a small configurational 
distance, but the opposite is not true. 

We have chosen silicon and magnesium clusters for this 
study since they are among the most extensively studied 
clusters and since we wanted to see whether clusters made 
out of insulating and metallic materials behave in the 
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The two structure were rotated and shifted in such a 
way as to minimize D. In addition atomic numbers were 
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III. RESULTS 

A. Silicon Clusters 

For silicon system we did our calculation for small clus- 
ters containing 3-19 atoms and for Si%2 as an represen- 
tative of medium size clusters. For very small clusters 
there exist only a few local minima structures and they 
are therefore usually well separated in energy. As the 
number of atoms in the cluster grows, the number of 
meta-stable structures increases exponentially. The con- 
cept of a global minimum is already rather ill-defined 
for silicon clusters containing more than some 7 atoms. 
They have many quite distinct structures that are very 
close in energy to the global minimum structure [ll| . As 
a consequence more than one structure can be populated 
even at room temperature. A second consequence of this 
is that different density functionals can give a different 
energetic ordering of the various minima [12] and even 
with the most accurate Quantum Monte Carlo calcula- 
tions it is difficult to obtain the resolution necessary to 
predict the correct energetic ordering [ll|. In this study 
we are not claiming to identify the correct ground state 
structures of the studied silicon clusters, but instead we 
want to show general trends. Therefore we use standard 
density functional theory instead of the extremely ex- 
pensive Quantum Monte Carlo method. Considering the 
fact that completely different structures can be extremely 
close in energy suggests strongly that a major perturba- 
tion such as the addition or removal of an electron can 
change the energetic ordering of the structures. Older 
studies have in contrast frequently just assumed that the 
ground state structures of neutral and positively charged 
clusters are the same. In some more recent investigations, 
few cases were identified where the neutral and positively 
charged cluster were not 'related'. In an investigation, 
where silicon clusters with less than 20 atoms were in- 
vestigated [13|, Sis, 5^12, ^13, Siis and Sin were found 
as the exceptions were the ground state geometries of the 
cation differ from the one of the neutrals. In another in- 
vestigation of silicon clusters with less than 10 atoms [lij , 
the ground state geometry of Sig and Si±o were found 
to be the "related" . Both studies are in contradiction to 
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FIG. 1: Global minima of charged and neutral Si n , for 
n=6,7,..19 and 32. Only for n=6,7,9 and 18 the global 

minima of charged and neutral are "related" .The 
configurational distance between each pair is given in A. 



our results which show that for silicon clusters with more 
than 7 atoms, the ground state structures of the neutrals 
and cations are not related with the only exception of 
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FIG. 2: Mapping chart for Si\$. The configurational 
distance between the the neutral and charged ground 
state configurations is very small (0.04 A) and ionized 
ground state does relax into the neutral ground state 
when an electron is added. However the neutral ground 
state does not relax into the ionized ground state and 
therefore the structures are not 'related' according to 
our definition. This behavior is rather exceptional and 
was only found for Si\o ,6^12 ,^#25 and Mg$Q. For all 

the other unrelated structures neither the ionized 
ground state relaxes into the neutral ground state nor 
the neutral into the ionized one. 



Si 14 



Neutral 



Charged 




FIG. 3: Mapping chart for Si\±. The ground state 
structures are not related and are quite different (Fig: [J). 



Szq and Si is and are as a matter of fact quite different ( 
Fig. [1]). In another study of medium sized clusters [HI 
it was also found that in most cases the structures of the 
neutrals and cations are the same. Out of the medium 
size clusters we have only examined the 32 atom clus- 
ter for which we however also find different ground state 
structures. 

Fig: [2] and Fig: [3] shows the mapping chart which 
gives detailed information about the relaxation proper- 
ties upon addition and removal of an electron . We dis- 
tinguish between reversible and irreversible mappings be- 
tween pairs of local minima. The energies of all the struc- 
tures are measured with respect to the ground state en- 
ergy of the neutral system. Solid double arrow connect- 
ing lines denote reversible mappings and dashed single 
arrow connecting line irreversible mappings. The space 
group is given in the rectangular boxes and the numbers 
close to the the connecting lines give the configurational 



distance of the two configurations. A reversible mapping 
connects two structures which are related according to 
our definition. In an irreversible mapping, the cluster re- 
laxes from the i-th to the j-th local minimum when an 
electron is removed or added, but it relaxes to a struc- 
ture which is different form the i-th when the electron 
is given back or taken away again. Fig: [2] and Fig: [3] 
shows that both kinds of mappings are encountered fre- 
quently. The minima of the neutral and cation are related 
according to our aforementioned definition only if a re- 
versible mapping connects the two global minima. This 
case was never encountered for clusters of more than 7 
atoms except for Sig and Siis and the global minimum 
structures for the neutrals and cations are thus different 
except for Si n n=3 to 7 , 9 and 18 in this size range . 
The numerical values along the relaxation arrows in the 
mapping diagrams indicate the configurational distances 
in the relaxation processes. These distances are typically 
of the order of 0.03 A, and thus show that the distor- 
tion during the relaxation is rather small. The symmetry 
group is also conserved in most cases. The fact that the 
geometries change so little upon removal or addition of 
an electron might have contributed to the wrong believe 
that the ground state of the neutral and cation are more 
or less identical. Nevertheless these small displacements 
are frequently sufficient to bring the system in another 
catchment basin. 

The energetic ordering for neutral and ionized clus- 
ter configurations would be identical if the ionization 
energy or electron affinities (including the energy that 
comes from the small relaxation upon removal or ad- 
dition of an electron) would be constant, i.e. indepen- 
dent of the shape of the various meta-stable configu- 
rations. The essential point is however that ionization 
energies and electron affinities are about two orders of 
magnitude larger than the energy differences between the 
ground state structure and the next meta-stable low en- 
ergy structures. Relatively small differences in the ion- 
ization energies and electron affinities between the differ- 
ent configurations can therefore lead to a reversal of the 
energetic ordering of the local minima. The energy differ- 
ences between the ground state and the first meta-stable 
configuration is of the order of few ksT at room tem- 
perature and the energy differences between the higher 
meta-stable configurations are even smaller. 

We find small configurational distance values not only 
for the structural changes induced by the addition or re- 
moval of an electron but also between different local min- 
ima of the neutral and ionized clusters. The configura- 
tional distance between the first and second meta-stable 
configuration of the Si 14 cluster is for instance only 0.15 
A. Nevertheless the two local minima are separated by a 
barriers of about 1.2 eV. In these disordered structure a 
broad distribution of barrier heights is to be expected [l6| 
and we find indeed also low barriers. The configura- 
tional distance between the ground state of the charged 
Si 10 cluster and its first meta-stable configuration is for 
instance also 0.15 A. But the barrier between the two 



structures are much smaller namely 0.22 eV and 0.08 eV 
respectively. Such small barrier heights are well below 
the accuracy level of density functional methods and it 
can hence not be excluded that higher level calculations 
such as coupled cluster or Quantum Monte Carlo cal- 
culation would give a different potential energy surface. 
Our previous experience [17j shows however that barrier 
height are quite well reproduced by density functional 
theory if no bonds are broken during the transformation 
from one structure to the other. 



B. Magnesium clusters 

For Mg n we have systematically studied all small and 
medium size clusters with n=6 to 30 atoms as well as 
M^56.The Global minima are shown in Fig. |4] and Fig. [5] 
For these cluster sizes the electronic HOMO-LUMO gap 
does not yet tend to zero, but is around 0.1 eV. So no 
pronounced metallic behavior is present. The ionization 
energies are also comparable to the case of the silicon 
clusters. The ionization energy is on average 5 eV for 
the magnesium clusters and 7eV for the silicon clusters. 
The only notable difference we found between the silicon 
and magnesium cluster is number of meta-stable states, 
which is much larger for silicon clusters. Since all energy 
differences are however also smaller for Mg than for Si 
,the average configurational distance between different 
meta-stable configurations is again similar in both cases. 
Hence Mg clusters have the same overall behavior as the 
Si clusters, i.e in general the neutral and ionized ground 
states are not related. 

In the studied size range we find that the global 
minima of neutral and cation clusters are related for 
n=7,8, 12, 15, 17,18, 19,24,26,27,30 and 32. For a bigger 
system ,Mg 5 6 we also found the global minima to be 
different for charged and neutral system. So in total the 
ground state structures are related in 12 cases and unre- 
lated in 21 cases. Fig. [6] and Fig. [71 exemplifies the same 
kind of mapping for MgiQ and M#24 between charged 
and neutral system as we already showed for silicon sys- 
tems. These mapping charts look very similar to that of 
silicon systems, i.e. the energetic ordering changes when 
the system goes from the neutral to the charged state. 
Although for M#24 the neutral and charged global min- 
ima are 'related' , from the mapping chart (Fig. [7]) we can 
see the sign of energetic ordering changes in the system 
while going from neutral to charged state. The numer- 
ical values along the relaxation arrows in the mapping 
diagrams indicate the configurational distances in the re- 
laxation processes. These distances are typically of the 
order of 0.02 A, unlike Silicon systems where this value 
is 0.03 A, and thus show that the distortion during the 
relaxation is smaller than that of silicon systems. 

Our results are again overall in disagreement with the 
majority of previous publications. In one of the earli- 
est publication on this topic, where clusters with up to 
6 atoms were studied, identical ground state structures 
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FIG. 4: Global minima of charged and neutral Mg n , for 
n=6-19. Only for n=7,8,12,15,17,18 and 19 the global 
minima of charged and neutral are "related" . The 
configurational distance between each pair is given in A. 
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FIG. 5: Global minima of charged and neutral Mg ni 
for n=20-30,32,56. Only for n=24,26,27,30 and 32 the 
global minima of charged and neutral are "related" .The 
configurational distance between each pair is given in A 



were found for Mg§ and Mgj [l8j |. In a study of Mg clus- 
ter with up to 21 atoms, it was found that only for Mg% 
and Mg4 the ground states are different In another 
somewhat more extensive study in the range between 2 
and 22 atoms [20|, it was found that in addition also 
M#6, Mg-?, Mg$, Mgn, Mg\ 2 and Mg\s have different 
ground states. 

We have also recalculated the energetic ordering of the 



minima of several magnesium clusters with the PBE func- 
tional [2l|. In all these cases the ordering was identical 
to the ordering with the LSD functional. This is in con- 
trast to the silicon clusters where the energetic ordering 
depends on the functional being used. This suggests that 
the density functional results for the magnesium clusters 
are very reliable. 

For the magnesium clusters the average configurational 
distance between the various local minima is typically in 
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FIG. 6: Mapping chart for MgiQ. The ground state of 
the neutral cluster is mapped to a rather high local 
minimum of the charged cluster. 
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FIG. 7: Mapping chart for M#24- For this system the 

ground states are related. The higher energy 
meta-stable states are however even for such a system 
typically not 'related'. 



the range between 0.1 A, and 0.2 A, and thus larger than 



the average configurational distance of the relaxation in- 
duced the the removal or addition of electrons. Since the 
magnesium cluster are also disordered we find, as in the 
case of the silicon clusters, a broad distribution of barrier 
heights. We calculated randomly 12 barrier heights of the 
neutral MgiQ and M#24 cluster and we found values in 
between 0.05 and 0.8 eV. 



IV. CONCLUSION 

Using an exhaustive sampling of the low energy config- 
urations based on the minima hopping method we show 
for silicon and magnesium clusters that the ground states 
of neutral and ionized clusters are in general not related 
and are in many cases quite different. This comes from 
the fact that for medium and large clusters there are in 
general numerous meta-stable structures which are ener- 
getically very close to the ground state. The differences 
in ionization energies and electron affinities for different 
structures are much larger than this energy difference be- 
tween structures. These facts have to be taken into ac- 
count in the interpretation of experiments with ionized 
clusters. 

There is no reason to believe that clusters made out 
of other elements behave differently. Based on our argu- 
ments one can only expect that for certain magic clus- 
ter sizes, for which ground state structures exist that are 
considerably lower in energy than other competing meta- 
stable structures, the ground state does not change upon 
removal or addition of an electron. Such an example is 
for instance the Cqq fullerene. 
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